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Abstract 

We examine the properties of open quantum systems with respect to 
their time evolution in different regimes, Markovian and non-Markovian. 
We analyze their behaviour with respect to their coherent or decoherent 
time evolution by means of different models and try to gain some insight 
into the possible correlations between Markovianity and coherence. 
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1 Introduction 

The properties of open quantum systems has a long history which starts with 
the measure problem already at the eve of the existence of quantum mechanics 
[DIO SO l[7l 11 HD]- The question of the measure of a quantum observ¬ 
able with a device which works on a classical level is related to the correlated 
concepts of decoherence and irreversibility. Although much progress has been 
made there remain many open questions leading to critics and controversies on 
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the subject m- An answer to these questions is of a particularly strong interest 
nowadays to the community which develops an intense activity in the framework 
of new fields like information theory and quantum computing. 

The central property of an open quantum system is the characteristic be¬ 
haviour of its time evolution. There are two aspects concerning this point. The 
first aspect is linked to the length of the time interval over which the system 
keeps the memory of its interaction with its environment. If this interval is 
short compared to the proper characteristic evolution time of the system the 
memory of the interaction is quickly lost, the regime is called Markovian. If 
this is not the case memory effects induced by a finite time response of the 
environment are present, so called backflow effects are generated, the process is 
called non-Markovian. 

The second aspect which is closely related to the measure problem concerns 
the evolution of the components of the density operator in a fixed basis of states. 
Given a specific basis of states of the system the evolution can either lead to so 
called decoherence if the non diagonal matrix elements of this operator decrease 
to zero, coherent if this is not the case. The question is whether and how the 
Markovian or non-Markovian property is correlated to a coherent or decoherent 
behaviour of the open system. 

Time behaviour and structure of the system-environment complex are nec¬ 
essarily intimately linked m- In fact Markovianity (non-Markovianity) and 
coherence (decoherence) are necessarily generated by the physical properties 
of the system, i.e. its environment and the interaction which couples them 
dynamically. This central point appears in the background of many different 
studies but has never been, to our knowledge, systematically analyzed as such. 
In the present work we aim to consider this point by means of models which 
allow to show how the structure of the different components of the total system 
determine the time evolution of the open system. 

The structure of the paper is the following. In section 2 we introduce a 
general formal expression of the density operator. In section 3 we develop two 
Markovian systems of different types which show different time evolutions. Sec¬ 
tion 4 is devoted to a non-Markovian system and in section 5 we show explic¬ 
itly how the algebraic structures of the Hamiltonian of the environment and 
the coupling interaction between the system and the environment determine 
the behaviour of the non-diagonal elements of the density operator. Section 
6 shows explicitly how the spectral properties of the environment can fix the 
transition from a Markovian to a non-Markovian behaviour and its consequence 
with regard to the type of evolution, either coherent or decoherent. Section 7 
summarizes and analyzes the results. Details of calculations are developed in 
Appendices. 


2 The density operator 

We consider an open quantum system S characterized by a density operator 
ps(t) which evolves in time from to to t under the action of the evolution operator 
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ps(t) = f(t,to)ps(to) 


( 1 ) 


At the initial time to the system S is supposed to be decoupled from its 
environment E and characterized by the density operator 



( 2 ) 


and in E space the density operator p E reads 

/5b( to) = ^2, dai : a 2 |op) (^ 2 1 


(3) 


ai 2 


Here |*i), \i 2 ) and |ai), |« 2 ) are orthogonal states in S and E space respec¬ 
tively, CjjjCjj normalized amplitudes and d ali0l2 weights such that p 2 E {to) = 
PE{to )■ 

At time t > to the reduced density operator in S space is ps{t) = Tr E [p(t)\ 
where p[t) is the density operator of the total system S + E. It can be written 

as [13] 



(4) 


with 



(5) 


where the super matrix C reads 



ol i ,a 2,7 


with 


U(i lh ),{a 17 ) (Mo) = 0’i7l^(Mo)l*iai> 
^(*4252),(«a7)(*>*«) = <*2«2|t/ + (t, t 0 )|j 2 T> 


(7) 


The unitary evolution operator reads U(t,to) = e lH< - t where H is the 
total Hamiltonian in S + E space and the super matrix C obeys the condition 


lim*—s-io ,i 2 ),(ji J 2 ) (^5 ^0) — Ja¬ 

in the present formulation the system is described in terms of pure states. 
The results which will be derived below remain valid if the initial density op¬ 
erator at the initial time is composed of mixed states ps(to) = XX i 2 c u*2 

We shall now use the general expression given by Eq.(4) in order to work out 
the density operator and matrix elements of ps(t) in the framework of different 
systems in order to characterize their time evolution properties with respect to 
Markovianity and coherence. 


3 Two Markovian systems 

The most general Hamiltonian acting in S + E space can be decomposed as 
H = Hs + He + Hse where Hse generates the interaction between S and E. 

It has been shown that the Markovian or non-Markovian property of S can 
be characterized by the presence of a unique state in E space mm or, more 
generally, by the commutator C'h = [He, Hse] [El- The evolution of S will 
be Markovian if Ch = 0 and non-Markovian if Ch 7 ^ 0. The consequences of 
these structures are analyzed by means of the following models. 

3.1 Model 1: The environment space reduces to a single 
state 

The Hamiltonian of the total system H = Hs + He + Hse in this order reads 


H = oosJz + f3a + a + p (a + J- + aJ+) ( 8 ) 

where J and J z are the angular momentum operator and its projection on the 
Oz and J_ = J x — iJ y , J + = J x + iJ y . The operators a + and a create and 
annihilate bosonic quanta, [3 and 77 are real parameters. 

Using the Zassenhaus expansion mm of the propagators U(t,to) and 
U + (t,to) which appear in the expression of the density operator it is possible to 
develop the exponential in terms of the Hamiltonian operators, see Appendix 
A. The expression of H is developed into its components and use in made of 
the commutation relations obeyed by the components of the angular momentum 
operator J. The development consists of an infinite series of terms which cannot 
be worked out in practice. 

However the structure of the development can be obtained. This is possible 
because there is only one state |a) in E space. In this case the only contribu¬ 
tions to the matrix elements of U. U + are the exponential terms which contain 
a + a,J z and J_J+ as well as their combinations. 
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As a consequence the matrix elements ps{t,to) are also diagonal in S space 
since the angular momentum operators defined above lead to diagonal matrix 
elements in the natural basis of states | jm) which are eigenvectors of J 2 and 
J z . The expression reads 


o) = \cj\ 2 \u ja {t,t 0 )\ 2 (9) 

where Uj a {t,t o) = (ja\U{t,to)\ja). In this example it comes out that there are 
no non-diagonal elements in the expression of the density operator whatever the 
size of the system space S, hence no decoherence because of the structure of the 
environment. In the sequel to = 0. 

3.2 Model 2: A Markovian system 

We consider another Markovian system for which [He,Hse\ = 0 in order to 
test its property with respect to coherence. Here 


Hs = wsa + a 
He = uJ z 

Hse = v(a + + a) J 2 (10) 

Using the corresponding notations introduced above the density operator 
ps(t) reads 


Ps(t) = X) X! C\nin2)(n3ri4) W c ni c n2 |n 3 )(n 4 | (11) 

ni,ri2 ri3,n4 


where the |rq)\s are now the normalized eigenstates of tfg. The super matrix C 
contains now the trace over the states [|ji?7ii)] in E space which are eigenstates 
of J z and J 2 , J 2 \jimi) = ji(ji + 1)| jirrii). 

Since Hs is diagonal in this basis of states C now reads 


)(n37i4) ($) ^{n^jm\U{t)\n l jm){n 2 jm\U*{t)\n A jm)/j i (12) 

jm 

and C'(niri2)(n^n4) (0) 1- 

Using the decomposition of the evolution operator e~ ltH already used in 
model 1 the projection on E space leads to 

= e ~ itul sa + a e -ir)j(j+l)sin(ust)(a + +d)/u Se u> s [l-cos(r)j(j+l)i\(a + -a ) 

The expression of C , (n 1 n 2 )(ra 3 ra 4 )(0 can now t> e explicitly worked out 
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( 14 ) 


C'(n 1 n 2 ){n 3 n4. ) (*>=£ -^ri2 n4 (*)4 3 ni W/5' 4 


with 


/„3„ 1 (t) = (n3|e- itw » a+ °e 

e ^ a [l-cos(r)j(j + l)t)]a+/j(j+l)? 7 e 


->13(i+l) sin(w a t)a+ e —i^'O'+l) sin(u; s t)a 

-«»[l-coB(»7jC7+l)t)o/jO'+l)»7] |m)e*(*) 


and 


( 15 ) 


In 2 rii 
^uj s [1—cos 


(t) = (n 2 \e +it “° a a 

(w'(l+l)i)]a + /10'+l)’7 


g+i^lO'+l) sin(a;,i)o + e +i?7j(j+l) sin(u; s t)a 
e -w»[l-cos(»jjC7+l)i)]o/j(j+l)7)| n4 \ e «(t) 


The matrix elements of /5s(t) read 


( 16 ) 


Ps 3 ” 4 ^) = E ^i< J C (« 1 " 2)( « 3 -. 4) ( t ) ( 17 ) 

ni,ri2 

The expressions of I n3ni (t) and I n2 n 4 (t) are explicitly worked out in Ap¬ 
pendix C by means of the Zassenhaus series development m, see Appendix A. 
They are exact. It can be seen that the expressions are oscillating functions of 
t , hence neither the diagonal nor the non-diagonal elements will decrease to and 
stay at zero when time flows. The system behaves coherently. 


4 Model 3: A non-Markovian system 

Consider the case where the Hamiltonian H of the total system reads 


H = Hs + He + Hse 


( 18 ) 


with 


Hs = ujJ z 
H e = (3b + b 

Hse = v(b + + b) J 2 ( 19 ) 
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which corresponds to the case where the role of S and E defined in model 2 are 
exchanged, 6 + , b are boson operators, oj is the rotation frequency of the system, 
P the quantum of energy of the oscillator and p the strength parameter in the 
coupling interaction between S and E. 

Now Ch 7 ^ 0. Since J z and J 2 commute in the basis of states [\jrn}\ the 
matrix elements of H in S space read 


(jm\H\jm) = urn + (3b + b + rjj(j + 1 )(b + + b) (20) 

The expression of the density operator ps{t) at time t is then obtained 
by taking the trace over the environment states of the total Hamiltonian p(t) 
leading to 


p s {t) = Tr E p[t) (21) 

whose matrix elements read 


ji mi ,.72 m2 

Ps 


(t)=p , 0 imuhm2 m E (ji,j2,t) 


( 22 ) 


with 


Jim 1 ,j 2 m 2 

P 0 


(*) 


g[— m 2 )]t 

0 ’lj 2) 1/2 


(23) 


with ji = 2 ji + 1. The bosonic environment contribution can be put in the 
following form 


&E(juh,t) = / . —7 


1 E nX (j 1 ,t)E*„ (j 2 ,t) 


n n ' I- 
n—0 n' ,n 


[K!)( 


n "!)] 1 / 2 


(24) 


The results are exact. The Zassenhaus development formulated in Appendix A 
has been used in order to work out the expressions m- The expressions of the 
polynomials E n ^ n '{t) and A 1 *,, n {t) are developed in Appendix B. 


By simple inspection of the expressions in Appendix B it can be seen that 
the non-diagonal of p-g mi ’ j2m2 (t) may cross zero when t increases but oscillate 
and never reach and stay at zero whatever the length of the time interval which 
goes to infinity. Hence no decoherence will be observable in this case. 


5 Model 4: A measure model 

As a fourth example we consider now a model which was introduced by Cuc- 
coli in the framework of a study of the correlation between measurement and 
decoherence [14]. 
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5.1 The model 


Following the von Neumann scheme for quantum measurements we consider a 
system S in contact with an environment E considered as a measuring device. 
The wave function of the total system S + E is separable at time t = 0 and 
takes the form 


l*(0)> = \i) ® |7> (25) 

where |z), |y) are sets of states in S and E space respectively. The [|z)]’s con¬ 
stitute a so called preferred basis of states of the system S which has to stay 
unchanged in time during the measuring operation 


|i>l7>-Hi) IV) (26) 

where | 7 *) are states which are affected by the presence of S space. Similarly 
for any linear combination of states |z) 


^Ci|z)|7) ^ ^Ci|z)|7 l ) (27) 


Introducing explicitly time in this scheme the total wave function at time t 
is given by 


!*(*)> =E c ‘l i >l / r < (*)> (28) 

i 

Under the assumption that the measuring device does not change with time, 
i.e. the operator corresponding to the observable to be measured commutes 
with the total Hamiltonian H = Hs + He + Hse, the expression of the density 
operator ps{t) reads 


ps(t) = E N 2 I*X*I + E CiC^Xi'K/WlE*)) (29) 

* ,7 


with 


|7 < (t))=e-«^| 7 ) 


(30) 


Here 


H l E = (i\Hs + H e + HsE\i) 


( 31 ) 



is an operator in E space 


H E — ti^E + He + H EE (32) 

where the e,’s are the eigenvalues of Hs- It conies out that the diagonal con¬ 
tribution to ps(t) is time independent. Hence the time evolution of the density 
operator is governed by the second term which is non-diagonal. 

Working out the second term in the expression of ps{t) for fixed state I 7 ), 
|i) and |i ) leads to matrix elements of the form 


M" 7 , = ac*,\i)(i \e 




f) e it(Ey-E' 


V \'y\e lt ^' SE \ r y )( r y \e ( 33 ) 


where q and E 1 are the eigenvalues of Hs and H E - 


The central point of interest concerns the behaviour in time of the non¬ 
diagonal terms given by Eq.(29) since they signalize the coherent or decoherent 
behaviour of the system. 

We consider two different cases which correspond to two different types of 
interaction Hamiltonians H l SE . 


5.2 Case a): [He, H 1 se ] = 0, i ^ i’ 

This corresponds to a Markovian behaviour of S since the matrix elements M 
are diagonal in E space mmm- In this case 


= 


Cj C.;, 


\i){i 


,it(e 


<7|e ftfl 4 «|7)(7|c-^-[7) 


(34) 


i.e. the environment (measuring device) stays in the same state all along its 
evolution in time. The contributions of the sum over 7 consist of real phases 
which add up as a sum of oscillations. A priori this never leads to the annulation 
of the non-diagonal terms M “ 7 for any finite or infinite time interval even though 
the sum may reach zero at some interval of times. Hence the system remains 
coherent. 


5.3 Case b): [H e ,Hse\ = kI e , i ^ i 1 

Here I E is the identity operator in E space. Using the Zassenhaus decomposi¬ 
tion [TSj the non-diagonal contribution of the super matrix M reads 


M"' 7 , = Ci4\i) (1 ) < 7 | e W^ E | 7 ') ( 7 ' \ e -iiHb E | 7 ) e «U/2 


(35) 


For the physical case where k < 0 the non-diagonal matrix elements of ps(t) 
decrease non monotonically to zero with a time scale 1/k. 
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6 Model 5: Markovian limit of an open system 
and decoherence 

As a last example we introduce a model proposed in refs. mmm in or¬ 
der to show how decoherence can be generated under certain conditions in the 
Markovian limit of a non-Markovian process. 

6.1 The model 

The system S' is a many-body system of particles with eigenstate energies [e*, * = 
1, N] coupled to an environment E of non interacting bosons or fermions through 
an interaction Hse characterized by a spectral density J a ij where the index a 
designates the particles in the environment. The spectral function reads mm 


J aij = v ^iKkA u - ( 36 ) 

k 

where V a ki is the interaction between the system states and the environment 
states |a). By means of the coherent state path integral techniques [25] in 
the framework of the Feynman-Vernon influence functional formalism 001 it is 
possible to derive the master equation which governs the evolution of dps{t)/dt. 
This equation has the structure of a Lindblad equation [21] with two terms. The 
second term is governed by time dependent dissipation coefficients They 

are obtained from the Green’s functions Gij(r,t 0 ) = ([a+(r), dj(fo)]). 

In matrix form 


c(t, t 0 ) = -l/2[dG(r, t 0 )/dtG 1 (r,t 0 ) + h.c] (37) 

6.2 Solution in a specific case 

The N * N Green’s functions G(t, to) obeys the following equation 

-j-G(t, f 0 ) + ?'e s G(r,t 0 ) + [ dT , v(r,r')G(r',t 0 ) = 0 (38) 

dr Jt 0 

where e s is the diagonal N * N eigenvalue matrix of the states in S and the 
N * N propagators v(r, r ) are given by 

' v (t,t') = ^2 [ ^Ja(w)exp(-iw(r-r')) (39) 

In order to be able to work out the solution of eqs.(38) and (39) in terms of 
explicit algebraic expressions we consider a diagonal spectral function J(w). The 
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uj dependence of J leads to dissipation functions Cu (t) which are time dependent 
functions and the behaviour of ps{t) is non-Markovian. 

If however J = J 0 = J 0 1 > 0 the propagator v(r, r ) = Jo<5(t — r ) and the 
solution of eq.(38) reads 


G(r,t 0 ) = exp-i(e s - iJ 0 )(r - t 0 ) (40) 

Introducing this quantity on the r.h.s. of eq.(37) leads to 


c(t) = Jo (41) 

hence to a constant quantity. As a consequence the master equation which 
governs the evolution of ps{t) gets a Lindblad expression corresponding to a 
Markovian process. In the general case the density operator possesses non¬ 
diagonal elements which monotonously decay to zero with time. 


7 Summary and conclusions: what can one learn 
from these examples? 

In the present work we investigated the coherence properties of open quantum 
systems evolving in a finite Hilbert space S and coupled to an environment evolv¬ 
ing in a finite Hilbert space E. We considered Markovian and non-Markovian 
systems and characterized this property by means of the structural properties 
of the different components of the total Hamiltonian of the system and its en¬ 
vironment. The investigation was performed by means of five models which 
encompass different physical cases. 

• The first two models described Markovian systems. In the first one the 
environment space E contained a unique state [12]. It came out that the 
system S was constrained to evolve in time in the state in which it was 
created at the origin of the process. Hence in the initial chosen diagonal 
basis of states the density operator of the system did not possess non¬ 
diagonal matrix elements and, as a consequence, no decoherence process 
can be observed. 

The second model generalized the former case. It was chosen such that the 
Hamiltonian of the environment He and the interaction Hse between the 
system and the environment commute, [He, Hse] = 0, which character¬ 
izes a Markovian process [1 Bill2) . This specific property which corresponds 
to a less stringent constraint than the one imposed in the first case led to 
a density operator which possessed also non-diagonal elements. However 
the non-diagonal elements came out to be oscillating functions of time 
which again did not lead to a decay to zero with time. 

Such a process is different from the expected behaviour generically ob¬ 
tained from the solution of a Lindblad equation. 
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• The third system described a non- Markovian process which was here 
defined by the commutation relation [He,Hse\ 7 ^ 0. It came out that 
there is no decoherence to be observed, the matrix elements of the density 
operator oscillate in time. 

• The fourth model relied on the von Neumann approach a measure pro¬ 
cess pH] . There the density operator was decomposed into a diagonal and 
a non-diagonal part with respect to a preferrential basis of states in the 
system space S. Two cases were examined: 

- In the first case where [He, Hse] = 0 the non-diagonal contributions 
added up as a sum over the environment states. They may behave in¬ 
coherently in time and decrease to or close to zero but this decrease is 
not monotonous,the contribution may again increase over later intervals 
of time and repeat to infinity. 

- In the second case the non-diagonal matrix elements were governed by 
[He, Hse\ — I where I is the unity operator in the total space S + E. 
There one observes an exponential decrease in t 2 modulated by oscillating 
contributions. 

• The fifth model was chosen as a description of the transition from a Marko¬ 
vian to a non-Markovian process. We chose a model already developed in 
ref. [21] which describes the time evolution of an open system coupled to its 
environment through an interaction given in terms of a spectral function. 
The master equation which governs the evolution of the density operator 
was put in the form of a Lindblad equation [241 ;25, 2S]- In the general 
case the dissipation functions which enter the second term of the equation 
are time dependent m , the system follows a non-Markovian process and 
does not necessarily decohere. However by extending the spectral function 
to a continuum with an infinite number of states the equation went over 
to a genuine Lindblad equation for the density operator. 

From these results it is possible to draw the following conclusions: 

• Markovian processes characterized by the divisibility property leading to 
[He,Hse\ = 0 m E3 EU do not necessarily generate a decoherent 
behaviour of the open system. This result is in agreement with formal 
previous work, see mmm- 

• A non-Markovian evolution which is characterized by [He, Hse] 7 ^ 0 may 
or may not lead to a coherent or decoherent process. The nature of the 
evolution depends on the structure of the evolution operator e ltH where H 
is the total Hamiltonian of S, E and their interaction, i.e. the commutation 
properties of the Hamiltonians of the system, the environment and the 
interaction which acts between them. 

• The transition from a Markovian to a non-Markovian process can be in¬ 
duced by the nature of the energy spectrum of the environment E. In the 
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model which was worked out mmm the interaction was governed by 
a spectral function which has an infinite energy extension. In this limit 
the evolution gets Markovian, it is governed by a Lindblad equation which 
leads to a monotonous decoherence process whereas a non-Markovian pro¬ 
cess generates different regimes which oscillate or not to zero [22 ] . 


8 Appendix A: the Zassenhaus development 

If X = —i(t - t 0 )(H s + He) and Y = -i(t - t 0 )H SE 

e x+Y = e x ®e Y ® e -^ X ’ Y )/» ® e -^ X ’ Y W ® e~ c ^ x ’ Y )H '... 
where 


c 2 (X,Y) = [X,Y\ 

c 3 (X, Y) = 2[[X, Y], Y] + [[X, Y], X] 

c 4 (X, Y) = c 3 (X, Y) + 3[[[X, Y],Y],Y] + [[[X, Y\,X\,Y\ + [[X, F], [X, Y], etc. 

The series has an infinite number of term which can be generated iteratively 
in a straightforward way [Hji . If [X, Y] = 0 the truncation at the third term 
leads to the factorisation of the X and the Y contribution. If [X', F] = c where 
c is a c-number the expression corresponds to the well-known Baker-Campbell- 
Hausdorff formula. 

Remark: the Zassenhaus expansion has a finite range of convergence. The 
upper convergence limit of time is in principle given by [2D] 


t=l/2\n2/(\\H E \\ + \\H S E\\) (43) 

Here the different models which are developed are analytically integrable, hence 
the series can be formally summed up to infinity. 


9 Appendix B: The bosonic content of the den¬ 
sity operator 

The expressions of the bosonic contributions to the density matrix p J s irn (t) 
are given by 


En,n'Uut) = e~ m Y, E (-*) n+n *(-l) T,,+Ba - n4 

n>n,2,ri3>n2 ri3>n4,n/>n4 


n\n'\(n 3 \) 2 [a(t) n+n3 ~ 2n2 ][^(t) n3+nl ~ 2rii ] 
(n - n 2 )l(n 3 - n 4 )\(n 3 - n 2 )\(n' - n 4 )\ 




(44) 
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and 


K’, n &h) = e if>t E E *"” +n 3 (-l )" + " 2_n4 

n” >n 2 ,n 3 >n 2 ri 3 >n 4 ,n>n 4 

n”!n!( 7 z 3 !) 2 [a(t)"”+ ra 3 - 2 Ti 2 ][Cft) ra+ra3 ~ 2ra4 ] ^ (t) 
(n” - 712)1(713 - ^ 2)!( n 3 — ri 4 )!(n - n 4)! 

The different quantities which enter E nin >(t) are 


( 45 ) 


a(t) = 


7(ji)sin/3f 

P 


(46) 


C(*) = 


/3[1 -cos 7 (ji)t] 
i(ii) 


(47) 


7(ii) = vhiji + 1) 


(48) 


,t, ^ 1 r 7 2 (ii)sin 2 (/3t) , /3 2 (1-cos 7 (ji)t) 2 1 

*''*> = — 2 - w - + -7(1)- 1 


(49) 


and for 


y(t) = 


7C72) sin Pt 

P 


(50) 


m = 


(3[1 - cos 7 [j 2 )t] 
7 O 2 ) 


(51) 


7 O 2 ) = vh(j2 + 1 ) 


(52) 


^2 (t) 


1 r7 2 (J 2 ) sin 2 (/3t) /3 2 (1-cos 7 (j 2 )i) 2 

2 P 2 + 7 2 (j2) 


(53) 
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10 Appendix C 

The matrix elements I n3ni (t) and I n2 n 4 (t) read 


I nani (t) = e~ iUst {-l) ni+m2+mi {-i) n3+m3 d l {t) n3+m3 ~ 2rn2 

m2rri3m4 


d 2 {t) ni+m3 ~ 2mi 


(n 3 — m 2 )!(?ri 3 — m 2 )l(m 3 - 7714 ) 1(711 — 7714 )! 


( 54 ) 


and 


^n,2 U4 ( t ) ^ 

d 2 (t) 


-\-iuj s t 

P2P3P4 
™4+P3~P4 


\ (—1 ) P2 +Pi+ " 4 (+i) n2 +P3 d\(t) n2 +P3 ~ 2p2 

n 2 \n4 ] -{p3 ] -) 2 e^^ 

(n i -p 4 y.(p3 -P2V-(P3 — 4*4)K n 2 -P2)! 


(55) 


For fixed 711 , 712 , 713,714 the sums over (?n,;,p.j) indices run over all values such 
that all factorials are positive numbers or zero. 

The coefficients di(t),d 2 (t) read 


di{t) = r)j(j + 1) sin(w s f)/w s 


(56) 


d 2 (t) = w s [l - cos {r}j{j + l)t)]\/rfj{j + 1) (57) 


and 


V(t) = ~V 2 \j(j + l)] 2 sin 2 (io s t)/2 (lo 2 ) - w 2 [l - cos{rjj(j + l)t)f/2{j{j + 1 ))V (58) 
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